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1 INTRODUCTION 

The origin of ultra-high energy ( > 10 18 eV) cosmic rays 
(UHECR) remains a long standing mystery in high energy 
astrophysics. Few objects are plausible as potential sources. 
In the Milky Way's ~ fiG magnetic fields, for example, 
UHECR have gyroradii comparable to, or larger than, the 
size of the Galaxy. Moreover, the apparent isotropy in their 
arrival directions, at least up to 10 19 ' 7 eV, indicates that 



they must have propagated a sufficient distance to have un- 
dergone multiple scatterings. This clearly points to an extra- 
galactic origin of UHECR. 

Consideration of the possible sources for these parti- 
cles has resulted in si ze, kinetic en ergy, and magnetic field 
strength constraints (|Hillasl 11984). with it being possible 
to infer little more beyond this. Recent observations by 
the Pierre Auger Observatory (Auger) have indicated that 
UHECR with energies > 10 197 eV, appear to have statisti- 
cally significant (at the 3cr level) correlations with the local 
(< 75 Mpc) AGN distribution in the V eron-Cetty catalogue 
AGN population (|Abraham et al 1l2008l). No such correlation 
was found in HiRes events ( Abbasi et alj |2008l ) , although 
the location of the observatory in the northern hemisphere 
required analysis with a different AGN set from the cata- 
logue. Nevertheless, such a signal is particularly interesting 
since it resides in the energy range where the pro pagation 
of UHECR becomes limited by GZK interactions l|Greisenl 
ll966l ; IZatsepin fc Kuzminll 19671 ). to several hundreds of Mpc 
in distance; an effect which should make any such correlation 
easier to see. 

One of the properties of cosmic rays at ultra-high ener- 
gies is that, for <nG extragalactic fields, they are capable 



ABSTRACT 

We investigate the acceleration of particles by Alfven waves via the second-order 
Fermi process in the lobes of giant radio galaxies. Such sites are candidates for the 
accelerators of ultra- high energy cosmic rays (UHECR). We focus on the nearby FR I 
radio galaxy Centaurus A. This is motivated by the coincidence of its position with the 
arrival direction of several of the highest energy Auger events. The conditions necessary 
for consistency with the acceleration timescales predicted by quasi-linear theory are 
reviewed. Test particle calculations are performed in fields which guarantee electric 
fields with no component parallel to the local magnetic field. The results of quasi- linear 
theory are, to order of magnitude, found to be accurate at low turbulence levels for 
non-relativistic Alfven waves and at both low and high turbulence levels in the mildly 
relativistic case. We conclude that for pure stochastic acceleration via Alfven waves to 
be plausible as the generator of UHECR in Cen A, the baryon number density would 
need to be several orders of magnitude below currently held upper-limits. 

of propagating great distances (tens to hundreds of Mpc) in 
the intergalactic medium almost rectilinearly. One source of 
particular interest is Centaurus A, a FR I radio galaxy in 
the southern hemisphere which has been associated with a 
number of the Auger detections abov e 10 19 ' 7 eV: two with 
the centre itself (jAbraham et al.ll2008h and additi onal others 
with the extended lobes ( Moskalenko et al1 l2008) . However, 
whether Cen A is merely deflecting particles towards our line 
of sight, or is itse lf the site of acceleration is a topic of cur- 
rent interest (e.g. iKotera fc Lemoi nc 2008; Gorbunov et al.l 
l200St) . 



Within the AGN class of objects, several different ac- 
celeration sites present themselves as candidate locations 
for the acceleration of UHECR. Among these are the cen- 
tral AGN, the jet driving the lob e, and the large radi o 
lobes themselves (for a review see iBegelman et al.lll994l ). 
Furthermore, several acceleration mechanisms may be act- 
ing in parallel, with different processes dom i natin g depend- 
ing on the local conditions (Ric ger et al. I [20071 ). In par- 
ticular, while first-order Fermi processes may be occur- 
ring near shocks, shear or inductive acceleration can op- 
erate in the jet llOstrowskil 1 19981 ; iRieger fc Duffy! |2004| ; 
iLvutikov fc Ouvedl |2007F T~ and stochastic acceleration is 
likely to occur within the lobes. 

In this paper we consider the case for stochastic accel- 
eration being responsible for the production of UHECR in 
giant radio galaxies. The lobes of these radio galaxies are 
most likely highly turbulent environments, where the frac- 
tional energy density in the magnetic turbulence, SB 2 /Bq, 
is not necessarily a small parameter. Due t o its size and esti- 
mated field strength, it has been argued bv lHardcastle et al.l 
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(2008) that stochastic acceleration may be a possible mech- 
anism for accelerating UHECR above 10 19 eV in Cen A. It 
was assumed in this work that the magnetic energy density, 
Btot/^TT, was within 10% of the rest mass energy density of 
the thermal protons in the lobe 0.5n p m p c 2 . 

Earlier numerical studies of stochastic acceleration were 
mostly concerned with methods for solving t he particle 
transport equation (e.g IPark fc Petrosianl 1 19961 , and refer- 
ences therein), or the use of stochastic differential equations 
l|Achterberg fc Krullslll992l ). Since tracking large ensembles 
of test particles is computationally expensive, previous in- 
vestigations considered either the short time dependencies 
of diffusion characteristics on turbu lence level in a narrow 
dynamic range of fcmax/fcm in = 10 3 (jMichalek fc Ostrowskil 
1996; iMichalek et al.l Il99^ ). or rapid accel e ration in fully 
turbulent strong fields (e.g lArzner et al.l l2006h . Atten- 
tion has also been directed towards the influence of co- 
herent acceleration processes invoked by means of an 
electr ic field component parallel to the lo cal magnetic 
field (|Dmitruk et al.l 120031 ; lArzner et al.ll2006l ). While par- 
allel fields may be legitimately admitt ed via resistive pro - 
cesses in the MHD fluid approximation. lArzner et alj (2006) 
have pointed out that if an electric field is generated via 
the Lorentz transformation of time dependent Fourier modes 
representing magnetic fluctuations, an underestimate of the 
acceleration times may result as a consequence of finite local 
parallel electric fields. 

An attempt at following the turbulent acceleration of 
protons in radio lobes ove r many decades in energy has 
recently been presented in iFraschetti fc Melial l|2008l ). al- 
though the acceleration timescales resulting from their nu- 
merical scheme appear to disagree with quasilinear theory 
predictions by several orders of magnitude. Furthermore, 
their results, taken at face value, seem to be in conflict 
with the Hillas criterion, which limits the maximum par- 
ticle energy i? ma x to which a given source can accelerate, 
£ msx = Ze^ABcR, where Z is the charge number, (3a is the 
Alfven velocity in units of c, B is the magnetic field strength, 
and R is the size of the source. 

The goal of this paper is to use quasilinear theory to 
place constraints on the parameters required to accelerate 
protons or ions to 10 18 eV and above, assuming a homoge- 
neous turbulent radio lobe. Numerical simulations are per- 
formed to check the accuracy of quasi-linear theory results 
when SB 2 /B 2 and Pa are not small parameters. The out- 
line of the paper is as follows. In Section [5] we review the 
basic results of stochastic acceleration using quasilinear the- 
ory, and discuss limitations on the acceleration timescales 
for UHECR. Section [3] focuses directly on the acceleration 
of protons and electrons in the giant outer radio lobes of 
Cen A. In Section|4]we discuss the numerical model, demon- 
strating to what degree the results of quasilinear theory hold 
in highly turbulent regions with non-linear magnetic turbu- 
lence or mildly relativistic Alfven velocities. Our conclusions 
are presented in Section [S] 



2 PARTICLE ACCELERATION IN 
TURBULENT FIELDS 

The motion of a charged particle, with charge q, in electric 
and magnetic fields, E and B respectively, is determined by 



the Lorentz force 

^=q(E + /3xB), (1) 

where (3 is the particle velocity normalized to the speed of 
light. Any change in the energy of a particle originates from 
the work done by the i5-field. In most astrophysical environ- 
ments, the local highly conducting plasma shorts out any 
large-scale electric field in the local fluid frame. However, 
for particles which cross back and forth across shocks there 
is no global frame in which all electric fields vanish. This 
is the underlying principle of diffusive-shock and shock-drift 
acceleration. In turbulent regions, it is the small-scale elec- 
tric fields associated with the plasma waves that accelerate 

particles. 

It was shown by iTverskoil (|l967T ) that the evolution of 
the energetic particle distribution can be described by the 
momentum diffusion equation 



df(p,t) 
dt 
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which is valid for scattering times-scales much shorter than 
the crossing time of the physical system. 

The form of the diffusion coefficient due to interac- 
tions with hydromagnetic waves h as been discussed in de- 
tail in numerous publications (e.g. iKulsrud fc Ferrari|[l97ll ; 
ISchlickeiserll2002l ). We will focus on acceleration in a region 
of Alfvenic turbulence with a one-dimensional power spec- 
trum W{k) oc k- q , where SB 2 /8ty = J ft femax W(k)dk, with 
fc = 2n/\. The wavenumbers fc m in and fcmax correspond to 
to the longest (A max ) and shortest (A m i n ) wavelengths in the 
system respectively. The spectral index of the power spec- 
trum is denoted by q. 

Neglecting numerical factors, and using the wave spec- 
tr um described abov e, the momentum diffusion coefficient 
is (|Schhckeiserlll989D 







"- 1 P 2 c 2 q 
oc p . 



(3) 



where r g = pc/eBo is the gyroradius of a particle, and 
/3a is the Alfven velocity normalized to the speed of light. 
This result is valid for particles with gyroradii smaller than 
the correlation length of the field. The diffusion of particles 
with larger gyroradii are essentially independent of momen- 
tum since such particles interact with the entire spectrum 
(Tsyt OV1C hl l 19661 ). This results in an increase in acceleration 
times, a fact which we demonstrate numerically in Section[4] 
The associated mean free path, L, is 



Bl 

SB 2 V A ma: 



(4) 



Equation [2] contains a systematic acceleration timescale 
given by 



tacc — 



Dip) 



Pa — «. 



2-8 



(•>) 



For the rest of this section we will consider exclusively a 
one-dimensional Kolmogorov spectrum of Alfven waves with 
q = 5/3, and a turbulence level SB 2 /B 2 = 1.0. While Kol- 
mogorov turbulent spectra are typical for such environments 
the actual value of the turbulence level is more difficult to 
determine, particularly in radio lobes where the magnetic 
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Reynolds numbers are large. However, the results should 
be approximately accurate for gyroradii much less than the 
correlation length of the field since we calculate the Alfven 
speed using the field's rms value. 

For gyroradii larger than correlation le ngth of the tur- 
bulence the acceleration rate falls off as p' ~ 2 (|Tsvtovichl 19661 ; 
iBerezinskii et al.l ll990 

Hence, we consider the highest en- 
ergy protons to be those which resonate with the largest 
wavelength in the system, r 9 /A max = 1- The largest wave 
is taken to be an order of magnitude less than the size of 
the lobes A max ~ 0.1 -Riobc- We note that the gyroradius 

for a proton of energy e is r g w 10 ( 10 if eV ) (rfc) k P c > 
such that for a radio lobe of size i?iobe = 100 kpc, and with 
B — 1 /iG, protons of energy 10 19 eV have gyroradii equal 
to A max . These parameter are in a similar range to those 
expected in the lobes of Cen A. We also note that protons 
with energy 10 20 eV (in the same magnetic field) have gyro- 
radii comparable to the size of the system, making the giant 
lobes of Cen A marginal candidates for the acceleration of 
UHECR to these energies, provided the particle has had suf- 
ficient time to reach these energies. According to equation 

the corresponding acceleration time for such protons is 

^acc(£max) ^ A m ax 

Using Faraday rotation measurements upper limits can 
be estimated for the density in radi o lobes. These are typ i- 
cally of the order n p = 10 -4 cm -3 |Begelman et al.lll994f ). 
although the proton density can in fact be much lower. 
Adopting this upper value for the density in the lobe, the 
field and matter energy densities are roughly 0.025 eV cm -3 
and 5 x 10 4 eV cm respectively, indicating an Alfven speed 
of /3 A = (U B /U P ) 1/2 = 7 x 10" 4 where U B = B 2 /8ty and 
Up = 0.5n p m p c 2 . Using these values, the acceleration time 
of a proton with energy e < e m ax = 10 19 eV is 

which is longer than the Hubble time. Clearly, if ultra-high 
energy protons are stochastically accelerated in radio lobes, 
the Alfven speed must be considerably larger, ie. lower den- 
sity and/or higher magnetic fields. In fact, as we will demon- 
strate, for Cen A we require mildly relativistic Alfven speeds 
(5 A ~ 0.1. 

We have so far neglected particle interactions with mag- 
netosonic waves and acceleration due to transit-time damp- 
ing (TTD). It is well known that fast mode magnetosonic 
turbulence, when of comparable intensity to the Alfven com- 
ponent, can present a more effe ctive mechanism for particl e 
acceleration through TTD (e.g ISchlickeiser fc Millerlll998h . 
the acceleration time being shorter by a factor ln(/3^ 2 ). How- 
ever, this is still not sufficient to explain the acceleration of 
UHECR using the above canonical values for conditions in 
the lobes. Also, at Alfven velocities /3a > 0.05 the accelera- 
tion times become comparable. 

Furthermor e , as h as been pointed out bv lEichlerl (I1979T ) 
and lAchterberd (|l98lT ) , acceleration by fast-modes can be- 
come ineffective when the plasma beta /3 > m c /m p , where 

1 Actually, the transition from t acc oc p 2 ~ q to t aC c oc p 2 is 
smooth and starts to occur at gyroradii less than the correlation 
length r g < A c , as shown in Section[4] Thus the results that follow 
can be considered as rough upper limits. 



(3 = Pth/Us, with Pth the sum of partial pressures. If the 
positron fraction of the plasma becomes significant, and the 
corresponding proton density drops, the electrons will dom- 
inate the total pressure. If /3 satisfies the above constraint, 
most of the wave energy goes into near-thermal particles 
and the fast-modes are subject to rapid Landau damping. 
Assuming reasonable parameters for radio lobes we find that 
/3 ~ 1 and therefore fast modes may be effectively dissipated. 
However, this remains a topic of on-going in vestigation (e.g. 
lYan fc Laza rian 200 j; [Petrosian et al.ll2006h . 

It may be the case that TTD makes a significant con- 
tribution to particle acceleration in localized regions within 
the radio lobes. However, a large departure from canonical 
parameters is required in order to accelerate particles to the 
highest energies and, as discussed above, it is likely that 
scattering from Alfvenic turbulence is the dominant process 
for the highest energy particles. In the investigations pre- 
sented here we have therefore focused exclusively on Alfven 
waves and defer to further studies the inclusion of fast mode 
effects. 

In Section U we investigate the accuracy of the acceler- 
ation timescales for Alfvenic turbulence provided by quasi- 
linear theory for both non-relativistic and mildly relativis- 
tic Alfven velocities, and for large 8B 2 /Bq. We emphasise 
here that the mildly relativistic models considered are highly 
speculative, involving ion densities three orders of magni- 
tude below the upper limits derived from Faraday rotation 
measurements. 



3 STOCHASTIC ACCELERATION IN CEN A 

Cen A is the nearest radio galaxy and one of the bright- 
est extraga lactic radio sources in the sky (for a review see 
llsraell fl998l . and references therein). Concentrated within a 
circle of 20° around the nucleus of Cen A, are 10 of the 
27 reported events above 10 19 ' 7 eV. While two appear to 
come fro m the nucleus itself, many more coincide with its gi- 
ant lobes ( Moskalenko et al. 2008), which are approximately 
~ 600 kpc in extent, and projected on the sky have an an- 
gular size of 10°. These correlations may well be deceptive 
with their sources possibly origin ating in the Centaurus su- 
percluster l|Gorbunov et al.ll2008T ). or simply scattered into 
our line of sight by the ma gnetic fields in the galaxy or lobes 
IlKotera fc Lemoine|[2008l '). 

Nevertheless, being the closest ra dio-galaxy it has long 
been favoured as a source of UHECR dFarrar fc Piran 2000; 
llsola et ahlbooift . It was investigated in lHardcastle et al.l 
(2008) whether these particles could be stochastically ac- 
celerated in the lifetime of the lobes. However, high Alfven 
speeds were found to be necessary, in agreement with our 
discussion in Section [2] The actual value of the Alfven ve- 
locity is uncertain. Using minimum energy arguments, the 
magnetic fields are not expected to exceed 3/iG in the giant 
lobes of Cen A (| Alvarez et al.ll2000r ). 

Due to uncertainty in the baryon fraction of the particle 
content within the lobes, an accurate measurement of the 
density is a much more difficult task, with little more than 
upper limit estimates possible. While different upper limits 
have been given in different regions of the lobes based on the 
soft X-ray thermal emission, there is majority agreement at 
~ 10 -4 cm -3 (for a review of different limits see sect. 7.3 of 
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lHardcastle et al" . 2008). However, these are only upper limits 
and the actual densities may be considerably smaller. 

Another uncertainty lies in the correlation length of the 
turbulence. While larger values allow for increased maxi- 
mum proton energy, smaller values may favour rapid elec- 
tron acceleration, provided the turbulent spectrum can ex- 
tend to sufficiently short wavelengths to resonate with the 
electrons. 

The turbulence is generated on large scales by hydrody- 
namic instabilities. Given that the projected physical extent 
of each lobe is approximately 250 kpc x f 00 kpc, with struc- 
ture evident on all observed scales, a reasonable estimate of 
the correlation length is ~ 10 kpc. 

The small-scale cutoff in the turbulent energy spec- 
trum occurs at the scale l eq at which equipartition be- 
tween the kinetic and magnetic energies is reached. The 
correlation length may be in turn related to i cq via 
the magnetic Reynolds number _R m according to i cq w 
^m 1,/2 'turb l|Ruzmaikin et al . 1989). Assuming a characteris- 
tic bu lk velocity of 0.1 c, and a temperature of 10 7 K. ISpitzerl 
( 1962) provides an expression for transverse resistivity which 
yields a magnetic Reynolds number R m ~ 10 35 . Hence, 
^cq <C r g thermal, where r g thermal is gyroradius of thermal 
protons. 

With A max w 10 kpc and B « 3 /iG, a maximum energy 
of £max ~ 10 19 ' 7 eV is indicated for r 9 < A ma x, which is in the 
energy range relevant for the Auger events. The acceleration 
time for a particle of energy e ma x is face ~ 10 5 /3^ 2 yrs, 
where we have assumed a Kolmogorov spectrum q = 5/3 and 
equipartition between the energy density in the turbulent 
and ordered magnetic field. 

The age of the radio lobes in Cen A is estimated to 
be 10 7 ' 5 yrs l|Hardcastle et al]|2008h . which suggests Alfven 
speeds /3a > 0.1 are necessary to accelerate particles to 
such high energies in the lifetime of the lobes. Again, as- 
suming 5 = 3 ^G as an upper limit indicates densities 
< iq-7.5 cm - 3 are necessary to achieve such values of /3a- 

3.1 Electron acceleration 

As an interesting aside, one can calculate the acceleration 
time of electrons in the lobes of Cen A. We assume a turbu- 
lent spectrum with a uniform Kolmogorov power law extend- 
ing to the length-scales required for resonance with radio 
emitting electrons. However, we note that the uncertainty 
associated with such an assumption is significant. 

The combined radiative cooling time for electrons via 
synchrotron in a 3 fiG magnetic field (Ub = 0.2 eV cm -3 ) 
and via inverse Compton in the CMB radiation field 
(C/cmb ~ 0.3 eV cm -3 ) is 

4 2 2 

tcooi = E e /-ca T /3 e j [/tot 

« 0.83 ( * V 1 ( — J) 1 Myr (7) 

V10 12 eV/ V0.5 eV cm" 3 / y v ; 

where [/tot = Ub + C/cmb- The corresponding acceleration 
time for electrons is 

Ucc ~ 25 (l0- 4 cm-J (iffkfc") 



Thus cooling and acceleration timescales for an electron be- 
come comparable at a critical energy 

3 / \ — - 

ni ( np \~ T I ^ max 1 2 
Sc ~ llO-*cm- 3 J VlOkpcj 

x( g»— ,W U «* a V 3 TeV (9) 

V0.2eVcm- 3 y VO.5 eV cm' 3 / w 

Below this energy, the acceleration or re-acceleration of elec- 
trons may not be safely neglecte d. 

Recently, lHardcastle et all ((2008) have presented re- 
sults of the radio spectra from selected regions of Cen A 
in the northern and southern lobes. These authors fit mod- 
els to the spectra in the southern regions using an inverse 
Compton cooled power law spectrum, where the spectral age 
increases further from the nucleus (both spectral ages being 
~ 10 7 ' 5 yrs) . In the northern regions the spectra obtained are 
different but can still be fit with a continuous injection bro- 
ken power law model. The northern middle lobe is claimed 
to be responsible for the observed asymmetric behaviour. 
The spectral breaks are found, in all regions, to occur in the 
range 10 9 — 10 10 Hz (corresponding to photon energies of 
10" 6 - 10" 5 eV). In the ~ ^G fields that are assumed in the 
lobes, this corresponds to synchrotron radiating electrons 
with energies > 10 10 eV, which falls below the critical value 
e c , although could be consistent with a lower magnetic field 
value. 

Interestingly, our value for e c was calculated using the 
upper limit for the number density. Lower densities such as 
those necessary for the stochastic acceleration of UHECR, 
increase the value of e c by several orders of magnitude, sug- 
gesting that stochastic acceleration may play an important 
role in the shape of the electron spectrum at these ener- 
gies. More self-c onsis tent calculations, sim i lar to those of 
iKardashevI (|l962f ) and lStawarz fc Petro sian (2008]), are war- 
ranted to investigate this further. 

Stochastic acceleration may also play an important role 
in the suggested re-acceleration of electrons observed near 
the h ot-spots of other radio galaxy jets (e.g iMeisenheimer] 
Il996h . 

Finally, we reiterate that, since the acceleration of 
10 10 eV electrons in a 3 pG field requires a Kolmogorov 
magnetic spectrum to span more than 9 decades in length 
scale with a uniform power law, any softening of the turbu- 
lence spectrum at large k will reduce the effects of stochastic 
acceleration. Hence, the large Alfven speeds required to ac- 
celerate UHECR, as found in the previous section, are not 
necessarily in conflict with the radio spectra. 



4 NUMERICAL TEST-PARTICLE 
SIMULATIONS 

The quasi-linear results of ISchlickeiserl |l989) used in the 
previous section are calculated assuming purely resonant 
interactions with a continuous spectrum of non-relativistic 
Alfven waves, /3a <§C 1, with total energy density [/ tU rb «C 
Ub- Since our results have taken these energy densities to 
be comparable, it is important to investigate the accuracy 
of the resulting expressions for the acceleration time. 

Numerical investigations of stochastic acceleration in 
Alfvenic turbulence have previously been performed by 
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iMichalek fc Qstrowskil (|l996h and iMichalek et~ai1 (|l999h . 
Their results show good agreement with quasi-linear the- 
ory although the simulations were limited in dynamical 
range, focusing on a single energy in a narrow range of wave- 
lengths. As yet, there has been no in depth numerical study 
of the energy dependence of the acceleration times. 

We calculate the acceleration timescales of test particles 
in low frequency Alfvenic turbulence, using an a priori field 
configuration which allows for time-dependent waves. 

Two different approaches are adopted to determine 
these timescales. 

Direct simulation is carried out if acceleration is suffi- 
ciently rapid that the initial particle ensemble can be evolved 
from (7} = f 5 to f 9 within 4 x f 7 gyroperiods (measured 
at the initial energy). This is the case for the mildly rela- 
tivistic turbulence models we consider in this work. 

If the acceleration timescales are very long such that 
particle ensembles cannot be evolved over the range of in- 
terest in computationally viable simulations, we take snap- 
shots of mono-energetic particle populations over the range 
and rely on the slow dispersion of the injected particle delta 
function distribution to derive early-time diffusion coeffi- 
cients for Gaussian fits to the particle distribution. This is 
the case for the non-relativistic models. 

We remark here that our investigations are limited to 
protons within the sample range (7} = 10 5 to 10 9 . Since we 
wish to evaluate the plausibility of Alfven mediated stochas- 
tic acceleration as a mechanism for driving protons to ultra- 
high energies, viable acceleration timescales over this range 
are a minimum necessary condition for acceleration from 
still lower energies. 

We note that protons with energies less than 100 TeV 
will undergo stochastic acceleration provided they have suf- 
ficient energy to satisfy the resonance condition with short- 
wavelength Alfven waves. This condition is necessary to fa- 
cilitate the pitch angle scattering required for maintenance 
of a quasi-isotropic distribution. 

We note that the accelerated particles require some pre- 
acceleration process to lift them from the thermal pool. 
However, the consideration of pre-acceleration is beyond the 
scope of this paper. 



4.1 Constructing the fields 

We have extended the test-particle transport code LYRA 
l|Q'Sullivan et al~l [2009) to include time-dependent mag- 
netic fields and corresponding time-dependent electric fields. 

A straightforward Lorentz transformation of magnetic 
fluctuations will, in general, lead to non-zero values of SEB. 
For the low frequency waves we consider, this is unphysical, 
and can lead to an unintended secular accele ration of the 
parti cle parallel to the local magnetic field l|Arzner et al.l 

l20oeh . 

In this work, we shall consider two distinct field con- 
structions which explicitly maintain 5E ■ B = 0. The first 
is a one-dimensional field (in the sense that k is parallel or 
anti-parallel to Bo). In the second construction, k is isotrop- 
ically orientated in three-dimensions. 

In each case, linear Alfven turbulence is imposed upon a 
mean field, Bo, by superimposing a large number of linearly 
polarized modes, SBj (j = 1..JV), according to 



SB j = Aj(B () x kj)e : 



i(kj -r-Uji+^j) 



(10) 



where kj and 



>j are chosen randomly within the dimension- 
ality constraints imposed upon the construction of the field 
as described below. The turbulent component of the field 
then app roaches homoge neity and isotropy in the limit of 
large N (|Batchelorl [19821. The random selection of Fourier 
phases, (j>j, ensures the field has no large scale coherent struc- 
tures, and the amplitudes Aj are chosen to create a Kol- 
mogorov power spectrum. The waves are taken to be purely 
Alfvenic such that 



= Pa I cosS I 



(11) 



with 9 the angle between the wave- vector kj and Bo. Since 
the direction of kj is determined using two angles in spheri- 
cal polar coordinates, chosen at random, the absolute value 
of the cosine must be taken to prevent a directional bias in 
the fields. 



4- 1.1 One- dimensional field 

In order to make contact with quasi-linear theory, we con- 
sider a simpl e one-dimensional f i eld co n struction similar to 
that used by iKulsrud fc Pearcd \ 19691 ). ISchlickeiserl (|l989T) 
and others, whereby weak perturbations have wave-vectors 
purely parallel and anti-parallel to the mean field. Explicitly, 
each wavevector is determined by a random sign selection 



kj d = ±B . 



(12) 



Vanishing SE ■ B is enforced by assuming a uniform 
linear polarization state for all modes. The polarization di- 
rection, 5Bi say, is arbitrarily assigned perpendicularly to 
Bo and then fixed for all remaining modes according to 



SBj = SB! V j. 



Via a Lorentz transformation the associated electric field is 
SEj = -^.A f (B x SB j )e 1 ^- r -^ t+ ^' > . 



(13) 
i\d is 
(14) 



4-1.2 Three-dimensional field 

Numerical investigations on the effe cts of anisotropy in th e 
turbulent field have been reported bv lMichalek et aD l|l999h . 
by considering waves drawn randomly from within cones 
of differing opening angle about the mean field. They find 
that the diffusion coefficients vary by less than an order of 
magnitude. 

However, it has also been a rgued by 

Goldreich fc Sridhaij (|l995T ): IChandranl (|2000r i; 
Cho fc Lazarian (2002) that Alfvenic turbulence is 
anisotropic, with eddies elon gated along the (l ocal) mean 



magnetic field (ie. k± > fcii). Yan fc Lazarian! |2004) have 



pointed out that this anisotropy can result in inefficient 
gyroresonant acceleration due to the reduction in the total 
power in parallel modes. The full importance of this effect 
is unclear in strongly turbulent fields, as may be prevalent 
in radio lobes. 

For instance, there is observational evidence from polar- 
ization maps for cha otic fields in the southern lobe of Cen A 
jjunkes et all 1 19931) . Furthermore, lobes such as those of 
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Cen A with mor phologies character ized by extended filamen- 
tary structu res (Carilli et al.l ll989) are known to be weakly 
mag netized (|Clarkelll996H ~ in the sense that SB > Bo. Con- 
sequentially, hydrodynamic turbulence, as might arise from 
high velocity injection of plasma into the lobes, will result 
in a strongly chaotic magnetic field component. 

We now describe the construction of a magnetic field 
from modes with wavevectors directed isotropically in three- 
dimensions and a corresponding electric field such that 
SE ■ B = 0. 

Appealing to Ohm's law for an ideal plasma with infi- 
nite conductivity, we may write 



SE = -Su x B, 



(15) 



where <5E, Su, B are the total electric, bulk plasma velocity 
in units of c, and magnetic fields respectively in the observer 
frame. 

Linearizing Faraday's law, and using the fact that linear 
Alfven waves are non-compressional (ie. kj ■ Suj — 0), we 
have 



Suj = — sgn(cos6*)^-5Bj 
Bo 



(16) 



for the contribution to the velocity field due to each Alfven 
mode. Hence 



5u = -§i(B+-B-), 
Bo 



(17) 



where B + (B ) is the sum of all forward(backward) moving 
modes with respect to the mean field. From equation ll5l then 



SE - 



-Su x Bo 



2^(B^ 
Bo 



x B 



(18) 



Therefore, retaining the second order correction to <5E given 
by the second term on the right-hand side of equation [18] 
provides a total electric field satisfying <5E- B = 0. In the ab- 
sence of this term, only SE ■ Bo = and the electric field can 
still have a parallel component to the perturbing magnetic 
field, leading to rapid electrical acceleration of a charged 
particle propagating in the total B-field, B. 



4.2 Equations of motion 

Given the time dependent fields B(r) and 5E(r), the rela- 
tivistic equations of motion for a particle of charge number 
Z, velocity (3 = v/c, and rest mass m are 



d7/3 
dt 
dr 
dt = 



Ze 



(8E + f3x B) 



/3c, 



(19) 



(20) 



where e is the electronic charge. In a magnetic field of mean 
magnitude -Bo, the characteristic relativistic Larmor angular 
frequency is ui g = ZeBo/jmc where 7 = (1 — /3 2 ) -1//2 . The 
corresponding maximal Larmor radius is r g = u~ j3c. 



4.3 Simulations 

To test the results presented in Section [2] we consider a 
homogeneous region of infinite extent with a mean field of 
Bo = 3 fiG onto which is superimposed an isotropic tur- 
bulent component SB. The turbulence level of the field is 



Table 1. Model parameters 
Model Dimension (5B/B ) 2 n p (cm -3 ) j3 A 



Wl 


1 


0.1 


10- 4 


0.002 


W3 


3 


0.1 


10- 4 


0.002 


XI 


1 


1.0 


10- 4 


0.002 


X3 


3 


1.0 


10- 4 


0.002 


Yl 


1 


0.1 


10- 8 


0.2 


Y3 


3 


0.1 


10-8 


0.2 


Zl 


1 


1.0 


10-8 


0.2 


Z3 


3 


1.0 


10-8 


0.2 



determined via the parameter {yj (SB 2 ) / Bo) 2 = (SB / Bo) 2 ■ 
The dynamic range of the turbulent field extends from 
A m in = 10-8 kpc to A max = 1 kpc and is resolved via 513 
modes evenly spaced logarithmically. We have confirmed 
that our results are adequately converged at this mode den- 
sity and dynamic range. 

The Alfven speed in the simulations is determined by 
the proton density according to 



(3 A =~ 



B„ 



c -y/47rra p n p 



(21) 



where m p and n p are the proton mass and number density 
respectively. 

We consider 8 different models in this work as defined 
by the parameters in Table [T] We consider non-relativistic 
turbulence Pa = 0.002 for both low and high turbulence lev- 
els (SB/Bo) 2 — 0.1 and 1.0: models W and X respectively. 
We also consider mildly relativistic turbulence (3a ~ 0.2 for 
{SB/Bq) 2 = 0.1: model Y ; and for (SB/Bo) 2 = 1.0: model 
Z . In each of these cases both one-dimensional and three- 
dimensional Alfven turbulence is considered. 

For the non-relativistic models W and X , the accel- 
eration timescales are too long to capture the evolution of 
an injected particle distribution over a significant energy 
range with sufficient accuracy to resolve any statistical sig- 
nal from numerical error. To see this, we have carried out 
benchmark integrations using a fifth order adaptive Runge- 
Kutta method with the accuracy tolerance set to 10- 6 . 

For model X3 , it takes approximately one hour on 
an Intel Xeon 2.66 GHz CPU to achieve an integration over 
6 x 10- 4 T acc (7 = 10 9 ) for a proton at 7 = 10 9 , or 4 x 
10" 5 T acc (7 = 10 5 ) for a proton at 7 = 10 5 . 

Therefore, for the non-relativistic models W and X , 
we rely on snapshots of the acceleration rate at energies 
l°Sio7 = {5> 6, 7, 8, 9}. Furthermore, we find that inte- 
grations over T = {64, 32, 16, 8, 4} x 10 3 Uo~g respectively, 
where ojo g (7o) is the relativistic Larmor frequency at the 
corresponding initial energy 70, are sufficient for the mo- 
mentum distributions to complete the initial ballistic trans- 
port phase and become diffusive. While t <§C r acc , the mo- 
mentum distributions may be approximated via Gaussian 
profiles. This can be easily demonstrated for the special 
case of a flat wave spectru m (q = 1) for which a solutio n 
to equation [2] is provided in iMichalek fc OstrowskH l|l996h . 
We consider termination times for integration ranging from 
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10" 



10° 



a 10° 



10 4 




10" 



10° 



10° 



10' 



10° 



10 3 



energy for three-dimensional turbulence with Runge-Kutta toler- 
ance parameter e = 10 — 6 . r acc is derived from systematic drift of 
momentum distribution via (A-y)fAt. Simulation results for 765 
particle ensembles at 5 different initial energies are indicated by 
filled circles. A best fit power law, 9O87 ' 69 is shown as a solid 
line. The expected values from quasi-linear theory for each of the 
case studies in this work are also represented via 

C-yl/3. W-long 

dash - C = 10224; X - short dash - C = 1022; Y - dot - 1.022; Z 
- dot-dash - C = 0.1022. 



4 x lCP 4 T acc (7 = 10 5 ) for particles injected at 7 = 10 5 , to 
1 x 1CT 2 T acc (7 = 10 9 ) for particles injected at 7 = 10 9 . 

It has been confirmed that energy is conserved to high 
accuracy by the integrator by examining its behaviour over 
characteristic integration times in the special case /3a = 
where r acc — * 00. Under these conditions for three- 
dimensional turbulence with (SB /Bo) = 1.0, energy is found 
to be conserved to within 0.005%. Figure [1] illustrates the 
systematic fictitious acceleration timescales experienced by 
particles in this instance. Reference lines indicating quasi- 
linear theory predictions for the acceleration timescales for 
each of the cases considered here are also shown. It is evi- 
dent that, assuming the numerical errors do not grow dis- 
proportionately for /3a > 0, sufficient numerical accuracy is 
employed to resolve the required signal from the statistical 
models. 

For each of the W and X model snapshot energies, 1500 
particles are injected with random position and velocity ori- 
entation in a randomly generated field. No more than three 
particles are released in any given realization in an effort to 
reduce numerical noise. 

The mildly relativistic Y and Z models are less 
challenging numerically. This is because the acceleration 
timescales are much shorter and therefore the integration 
can be carried out over a relatively long time. We inject 
500 particles at 7 = 10 5 for each of the one-dimensional 
and three-dimensional models and integrate continuously for 
4 x 10 Wrjg(7 — 10 ). This is sufficient for acceleration to 



values of 7 in excess of 10 



4.4 Non-relativistic turbulence 

As previously remarked, for non-relativistic turbulence with 
/3a = 0.002, integration over a large energy range is not 
feasible computationally. Instead, ensembles of particles are 
injected at discrete energies over the range of interest and in- 




-0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 

A7/70 

Figure 2. Model Wl : Histograms of relative energy change, 
A7/70 (where A7 = 7—70), for ensembles of 1500 mono-energetic 
protons injected at t = with 7 = 70- Solid line indicates a 
best-fit Gaussian distribution: the mean drift {A7} and variance 
{(A7) 2 ) are given in Tablc[2] 
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Table 2. Model Wl : Best-fit Gaussian distribution parameters. 



70 


V((7- 


<7» 2 >/70 


(A7)/7o 


10 s 


0.0107 


± 0.0005 


0.000506 


± 0.000583 


10 6 


0.0163 


± 0.0007 


-0.000290 


± 0.000897 


10 7 


0.0239 


± 0.0011 


0.000031 


± 0.001297 


10 s 


0.0343 


± 0.0015 


0.000887 


± 0.001795 


10 9 


0.0314 


± 0.0025 


0.000065 


± 0.003111 



10° 



o 10' 

u 



10° 



10 J 



~n — 1 

QLT 
Wl Sim 
W3 Sim 
Wl Fit - 
W3 Fit 



10° 



10° 



10' 

<7> 



10° 



Figure 3. Models W : Acceleration T aC c as a function of en- 
ergy. T acc is derived from best fits of Gaussian profiles to the 
standard deviation of the momentum distributions via 7q/-D p 
where D p = (A7 2 )/2At. The expected value from quasi-linear 
theory, 102247 1 / 3 , is represented by the solid line. Note that for 
p = ?'g/A c > 1 the efficacy of the acceleration process is damped. 
The best fit to the data for p < 1 ((7} < 3.6 X 10 8 ) is given by 
45317 ' 36 (dashed line) for Wl , and 191387 ' 36 (dotted line) for 
W3 . 



stantaneous acceleration times are inferred. This is achieved 
by allowing the injected delta function to relax to a Gaussian 
distribution and fitting its variance. 

In Figure [2] snapshots of the particle distributions are 
illustrated for the Wl model. The fitting parameters and 
errors are provided in Table [2] for reference. Evidently, while 
the mean drift upward in energy is overwhelmed by the er- 
ror, the standard deviation is well captured to within ~ 5%. 
We have confirmed that the variance of the distribution 
does indeed vary linearly after a few scattering times, r s , 
ie. becomes diffusive. Furthermore, all simulated distribu- 
tions achieve this diffusive phase of evolution with r s < l c /c 
for p < 1 and r s > l c /c for the highest energy test cases. 

The momentum diffusion coefficient is approximated by 



D p = 



(A 7 2 ) 
2At 



(22) 



where we have used the fact that the mean drift in the dis- 
tribution is negligible over the integrated times ie. A7 = 
7 — (7) w 7 — 70. From this expression, the acceleration 
time, race, is then obtained using 



(23) 



10' 



.10° 



■ 10° 



10" 



QLT 
XI Sim 
X3 Sim 
XI Fit 
X3 Fit 

_d I u 



10° 



10° 



10' 

(7> 



10° 



10 3 



Figure 4. Models X : Acceleration r aC c as a function of en- 
ergy, race is derived from best fits of Gaussian profiles to the 
standard deviation of the momentum distributions via 7q/-D p 
where D p = (A7 2 )/2At. The expected value from quasi-linear 
theory, IO227 1 / 3 , is represented by the solid line. Note that for 
p = r g /A c > 1 the efficacy of the acceleration process is damped. 
The best fit to the data for p < 1 ((7) < 3.6 X 10 s ) is given by 



587 u 



(dashed line) for XI , and 957°' 50 (dotted line) for X3 



The results for the weakly turbulent models Wl and 



W3 are plotted in Figure [3] In both cases an index of 0.36 is 
found for r aC c using a best fit containing all points excluding 
that for the 10 9 ensemble. This is in close agreement with the 
theoretical value of 1/3 dictated by quasi-linear theory, al- 
though the index does appear to approach this value asymp- 
totically for lower rigidities. Notably the three-dimensional 
field presents acceleration times that are a factor of 4 larger 
than the one-dimensional field case. This difference occurs 
due to the cos 9 dependence of the electric field strength for 
oblique Alfven waves. The fact that the same normalization 
for the total turbulent power was used for both simulations, 
SB 2 /8ir naturally accounts for the longer acceleration times. 
Comparing the results from the one-dimensional simulations 
with the quasi-linear theory result, we note that the numer- 
ically determined acceleration time is faster by a factor of 
almost 2. This may be due to the fact that the particles 
in our simulations can resonate with higher harmonics not 
accounted for in quasi-linear theory, or simply that the tur- 
bulence level is not low enough. Nevertheless, the results in 
all cases are in good agreement to order of magnitude. 

It is also observed that for p>l (7^ 3.6 x 10 8 ) the ac- 
celeration time strongly diverges from the theoretical result 
as particles' gyromotions exceed the scale at which resonant 
waves are present in the field. 

Regarding the acceleration of UHECR, this emphasizes 
the fact that the maximum energy is quite typically governed 
by the correlation length of the field. 

For models XI and X3 , where the turbulence level 
is high, the agreement with quasi-linear theory is less com- 
pelling, as illustrated in Figure[4] Here, the indices found are 
0.46 and 0.50 for the one- and three-dimensional cases re- 
spectively, and the experimentally determined acceleration 
times in both cases are also an order of magnitude faster 
than predicted in quasi-linear theory. The divergence from 
quasi-linear theory values with increasing turbulence level 
is not surprising, as formally quasi-linear theory is valid in 
the limit of weak perturbations to the mean field. Again, the 
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t(A c /c) 

Figure 5. Models Y : Mean energy against time for ensembles 
of 510 particles injected at 70 = 10 5 . Reference line indicates 
p = 1. quasi-linear theory result is {7) = 0.03465t 3 . Yl numerical 
results are indicated by long-dashed line; Y3 numerical results 
are indicated by short-dashed line. The critical energy 3.6 X 10 8 
at which p = 1 is indicated by a dotted line. 



10" 




QLT 

Zl 

Z3 

p = 1 



_l ' 



10 J 



10* 



t(A c /c) 



Figure 6. Models Z : Mean energy against time for ensembles 
of 500 particles injected at 70 = 10 5 . Reference line indicates 
p = 1. quasi-linear theory result is (7) = 34.65t 3 . Zl numerical 
results are indicated by long-dashed line; Z3 numerical results 
are indicated by short-dashed line. The critical energy 3.6 X 10 s 
at which p = 1 is indicated by a dotted line. 



three-dimensional field shows acceleration times longer than 
the one-dimensional case by factor of ~ 4. 



4.5 Mildly relativistic turbulence 

For the mildly relativistic models Y and Z , with /3a = 0.2, 
n p has been set to 10~ 8 cm -3 , a value four orders of magni- 
tude below Faraday rotation upper limits. While this value 
may be extreme, we wish to demonstrate the conditions un- 
der which second-order Fermi acceleration via pure Alfvenic 
turbulence is viable in our idealized model. We note that 
such low densities require a large fraction of positrons in or- 
der to safely neglect the influence of fast modes, as discussed 
in Section [2] 

We find that for models Y and Z , the systematic ac- 
celeration is strong enough that we may simulate the accel- 
eration of a single ensemble of particles from a low injection 
energy, 70 = 10 5 , up to values of < 10 10 . Furthermore, fewer 
particles are required in order to derive a smooth signal. In 
practice this corresponds to integrations of 500 particles per 
model over 4 x 10 7 o;^ 1 (7o). 

In these cases, the one dimensional field models, Yl 
and Zl , exhibit more rapid acceleration than predicted 
by quasi-linear theory, with the low and high turbulence 
cases showing similarly curved profiles above the predicted 
(7) ~ t 3 power law. However, due to the relatively rapid 
acceleration in the plane perpendicular to the mean field, 
isotropy of the particle distribution cannot be maintained. 
As such, comparisons with the results of quasi-linear the- 
ory are misleading. 

The mean energy plots for the three-dimensional field 
models, Y3 and Z3 , have flatter profiles. The high tur- 
bulence model is in excellent agreement with quasi-linear 
theory and the lower turbulence level case has a somewhat 
slower growth rate than predicted by theory. 

As before, in all cases, once p>l, the acceleration is 
severely damped by lack of resonant scales in the field. 



5 CONCLUSIONS 

The stochastic acceleration of UHECR due to Alfvenic tur- 
bulence in the lobes of giant radio galaxies has been consid- 
ered. The work is motivated by the recently reported correla- 
tion between a number of the highest energy events observed 
by the Pierre Auger observatory and nearby AGN. Making 
simple estimates, based on the results of quasi- linear the- 
ory, we have investigated the maximum energies that can 
be achieved with the use of reasonable parameters in radio 
lobes. 

While the maximum possible energy that can be 
achieved in each lobe is primarily determined by the coher- 
ence length of the magnetic field in the lobes, the timescale 
to achieve this is largely determined by the Alfven speed. 
Assuming homogeneity of all quantities throughout the en- 
tire lobe, we have placed constraints on the value of the 
Alfven speed necessary to stochastically accelerate particles 
10 18 eV within the lifetime of the lobes. 

We have demonstrated that Alfven speeds in excess of 
0.1 c are required in Centaurus A to accelerate cosmic rays 
above 10 19 eV, although heavier nuclei could, in principle, 
achieve even higher energies. Also, should the composition 
of UHECRs with energies above 10 19 eV be dominated by 
heavier nucleii, the observed UHECRs could be accelerated 
in the lobes of Cen A with an Alfven speed an order of mag- 
nitude smaller. While the absolute value of the magnetic 
field strength in the lobes is quite well constrained by min- 
imum energy arguments to fall in the range ~ 0.1 — 3 fj,G, 
the densities referred to in the literature are merely upper 
limits, and indeed the proton number density could be con- 
siderably lower. 

Provided the turbulence spectrum steepens consider- 
ably at smaller wavelengths, such high Alfven velocities are 
not in conflict with the synchrotron radio emission from the 
non-thermal electrons. Conversely, should the Kolmogorov 
spectrum extend to sufficiently small wavelengths, the re- 
sulting rapid acceleration time of electrons might present a 
problem to the standard picture of an injected power law 
with a cooling break. 
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We have not allowed for transit-time damping (TTD) 
in the studies presented here since fast magnetosonic waves 
are not included in the turbulent field. While this simpli- 
fies the model and its interpretation, it may be justified by 
considering the efficacy of energy losses from these modes 
through heating of near-thermal electrons via Landau reso- 
nance effects. A comprehensive study of the linear theory of 
stochastic acceleration using a fully relativistic approach is, 
to our knowledge, still lacking, but a naive extension of the 
non-relativistic results to mildly relativistic phase velocities 
suggest that acceleration by Alfven waves dominates above 
/3a ~ 0.1. 

The numerical results presented in this work demon- 
strate purely stochastic acceleration of relativistic charged 
particles through interaction with linearly polarized Alfven 
waves. The ideal MHD constraint, (5E • B = 0, is triv- 
ially maintained for fields with perturbations directed par- 
allel to the mean field (fD) and by deriving an appro- 
priate electric field from inferred velocity perturbations 
via Ohm's law for an ideal plasma (3D). For both these 
one-dimensional and three-dimensional field configurations 
we consider low and high turbulence levels, (SB/Bo) 2 = 
0.1, 1.0, for non-relativistic and mildly relativistic turbu- 
lence, /3 A = 0.002, 0.2. 

For non-relativistic turbulence, we find that agreement 
with quasi-linear theory is good at low turbulence levels. The 
three-dimensional acceleration timescales are found to be in 
excess of quasi-linear theory by a factor of ~ 2. At higher 
turbulence levels, this agreement is less pronounced with the 
power law dependence of the acceleration time on the energy 
becoming steeper. Over the range considered in this work, 
the absolute acceleration times are also shorter than pre- 
dicted by quasi-linear theory with approximately the same 
relative difference between the one-dimensional and three- 
dimensional cases. 

Notably, above the critical energy corresponding to res- 
onance with the largest wave in the system, p = 1, the par- 
ticles interact with the the entire Alfven spectrum such that 
the diffusion becomes almost independent of energy. This 
dramatically reduces the efficacy of the acceleration process. 

For the simulations with mildly relativistic Alfven 
speeds, the mean energy of the particle ensembles shows 
reasonable agreement with quasi-linear theory up to p — 1 
for both low and high turbulence levels. The one-dimensional 
cases have more rapid acceleration than predicted by theory, 
but the distributions become anisotropic, and comparisons 
with quasi-linear theory are not applicable. 

Given that the results of the numerical simulations are 
in reasonable agreement with the results of quasi-linear the- 
ory, in the mildly relativistic regime, the conditions neces- 
sary for stochastic acceleration of UHECR in Cen A can 
be satisfied, provided the correlation length of the magnetic 
field is sufficiently large, and that a particle can remain in- 
side the lobe for a sufficiently long length of time. Paramount 
to this result, is the relative low density of baryons in the 
lobes, which must fall considerably below the upper limits 
inferred from observations. 
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